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ABSTRACT 

Expansion of non-spherical relativistic blast waves is considered in the Kompaneets 
(the thin shell) approximation. We find that the relativistic motion effectively "freezes 
out" the lateral dynamics of the shock front: only extremely strongly collimated shocks, 
with the opening angles A8 < 1/T 2 , show appreciable modification of profiles due to 
sideways expansion. For less collimated profiles the propagation is nearly ballistic; the 
sideways expansion of relativistic shock becomes important only when they become 
mildly relativistic. 



1. Introduction 



Dynamics and corresponding radiative signatures of non-spherical relativistic shocks remains 
an important unresolved issues in studies on Gamma Ray Bursts (GRBs). Since GRBs produce 
narrowly collimated outflows that evolve laterally, understanding the overall dynamics - both the- 
oretical and in terms of agreement between different numerical results - is imperative to the inter- 



pretation of the broadband observations of GRBs Rhoads (1999); Frail et al. (2001). 



Presently, there are two competing views on the lateral evolution of the relativistic outflows. 
Theoretically, it is typically argued that the lateral evolution of the flow proceeds with relativistic 
velocities (Piran 1999), (see also Wygoda et al. 2011). This view is contradicted by the results 
of numerical simulations that show very little lateral evolution in the relativistic regime |Cannizzo] 



et al. (2004); Zhang & MacFadyen (2009); Meliani & Keppens (2010); van Eerten & MacFadyen 



(2011). 



In this Letter we argue that this disagreement results from the incorrect theoretical assumptions 
about the lateral evolution of the flow. What is important for the interpretation of observations is 
the evolution of a curved shock. Previously the lateral evolution of the non-spherical shocks was 
incorrectly treated as a free lateral expansion into vacuum (e.g., Wygoda et al.pbll Eq. 5).. The 
assumption of the lateral expansion with the sound speed results in a "gramophone-type" profiles 
and exponential slowing down of the ejecta. This has drastic implications for the underlying light 



curves (eg Kumar & Panaitescu 2000 Panaitescu & Kumar 2003). In fact the dynamics of the 
non-spherical shocks is more subtle; the correct treatment, as we argue below, is consistent with 
slow lateral evolution seen in numerical simulations. 
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Evolution of strong non-spherical shocks is a well studies problems in fluid dynamics. The 
two fundamental works that have laid the foundation for non-spherical (two-dimensional) shocks, 



due to Kompaneets (1960) and to Laumbach & Probstein (1969), were originally designed to treat 



strong shock waves in the non-isotropic medium. These two complimentary methods have been 
extensively applied in astrophysics to treat supernova explosions ( Bisnovatyi-Kogan Sz Silich||1995 ) 
and non-isotropic winds (e.g., Icke|[l9 88), In the Kompaneets approximation the internal pressure 
of the gas is assumed to be constant. Then the Rankin-Hugonio conditions determine the normal 
velocity of the shock in the external inhomogeneous medium. A modification of the Kompaneets 



approximation - a thin or snowplow shell approximation - has also been used extensively (e.g., Wiita 



1978 Mac Low & McCray 1988 Bisnovatyi-Kogan et al. 1989). In a complimentary Laumbach 



Probstein approach (Laumbach & Probstein 1969) the streamlines of the shocked material are 



assumed to be radial, thus neglecting the lateral pressure forces. 

The relativistic generalization of the Kompaneets and the Laumbach-Probstein methods have 



been discussed by Shapiro (1979). Relativistic dynamics provide extra support for the thin shell 



method, since in the relativistic blast waves the shocked material is concentrated in even narrower 
region R/T 2 than in the non-relativistic Sedov solution. In addition, the limited causal connection 
(over the angle ~ 1 /T) provides a justification for the Laumbach-Probstein method on the angle 



scale comparable to 1 /T. As has been pointed out by Shapiro ( 1979 ) , the two methods - Kompaneets 



and Laumbach-Probstein - become very similar in the relativistic regime. This is due to the fact the 
in a relativistic quasi-spherical wave, the typical angle that a shock wave makes with the direction 
of the velocity is of the order a ~ 1 /T 2 . Thus the post shock pressure along the shock differs only by 
one part in T 2 , so that both approximations of constant post-shock pressure and radial post-shock 
motion become equivalent. In some sense the propagation of strongly relativistic non-spherical 
shocks becomes trivial: relativistic kinematic effect freeze out the lateral dynamics of the flow so 
that different parts of the flow behave virtually independently. 



2. Relativistic non-spherical shocks in the thin shell approximation 



In this section we re-derive the relativistic Kompaneets equation Kompaneets| ( |1960 ) ; Shapiro 



( |1979 ) allowing for the arbitrary velocity of the shock and arbitrary (angle-dependent) luminosity 



and/or external density. Consider a shock propagating with a three- velocity V at an angle a to its 
normal. There are three generic rest frames in the problem: laboratory frame K, a frame where the 
shock is normal to the flow K\ and a shock rest frame Kq. A frame K\ is related to the lab frame 
K by a Lorentz boost along y axis with a Lorentz factor T\\ = — V 2 sin 2 a. In K\ the velocity 

of the shock is V\ = T\\V cosa (along the x direction). Thus, r 2 = 1/(1 — V 2 ) = Y 2 cos 2 a + sin 2 a 
(the shock becomes non-relativistic when tt/2 — a ~ 1/r). In the frame Kq V\ and T\ are the 
velocity and the Lorentz factor of the unshocked medium. In the lab frame the x component of the 
shocked velocity V' x = V{/(T^(1 + V{V sma) generally has a completed form, but simple relations 
can be obtained in the strongly relativistic limit (see below). 
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We introduce next an acceleration parameter K (Kompaneets 1960, Icke 1988) as a Lorentz 
factor of the shock in the K\ frame 

K = T\ (1) 

The acceleration parameter can be expressed in terms of a ratio of a post shock pressure P' to the 
upfront density p. For relativistic strong shocks with the ratio of specific heats 7 = 4/3 

2 P' 

K = ^~ (2) 

3 p 

while for non- relativistic shock with 7 = 5/3 

Expressing the relevant quantities in terms of K we find 



|, ;K — sin 2 a y/K 



cos 2 a cos a 



K — 1 cos 2 a 

v/ if -sin 2 a ~ 2iT 



Vi = a/1- — ps 1- — 
V K 2K 



K — sin 2 a l 



K cos 2 a cos a 



, cosa / (2 + sina) cot 2 a/2\ 

x ~ 1 + sina 1 2if 1 ' l ' ;i 



where the approximations assume strongly relativistic motion. 
Finally, expressing V in terms of V\ we find 



V 



2_ V, 2 1-1/ K 



V 2 sin 2 a + cos 2 a 1 — sin 2 a /if 



cos- 5 a 
=2, 



if Vi <C 1 



S 2 ( 

1 - if Vi -> 1, arbitrary a 

+ ifa^O 

(5) 



Consider next a small section of the non-spherical shock at the spherical polar angle 9 propa- 
gating at an angle 

dlni? 

tana = ^— (6) 



to the radius vector (Fig. [T]). Then 



- 4 - 



(dRY 2 = ^-VK 

V 7 1 K(l+{9ehs.R) a ) 

Here K is a function of the shock position and angle K = K(9, R). Equation ([7]) is the sought rela- 
tivistic generalization of the Kompaneets equation (the thin shell modification of the Kompaneets 
equation, to be more precise). For example, for non-relativistic V = 1 — 1/K <^1 eq. ([7]) reduces 
to the familiar Kompaneets form 

— = vya + ^mi?) 2 ) (8) 

The other two simplifying cases of the Kompaneets equation include relativistic motion, Jf>l: 

dR\ 2 1 / 1 



dt J 1 K\l + (8 ]nR) 2 ) (9) 



and arbitrary quasi-spherical motion, a <C 1, if-arbitrary 

dR\ 2 ( . \ \ ( . 1 , „ 2 



m) V 1 "^l 1 + ^ lnfi) J (10) 

Again, the last equation readily gives the standard Kompaneets equation as K — > 1 and 1 — 1/ K — > 
V 2 . 



Examination of the eq. ([7]) confirms that for relativistic motion the angle a ~ 1/(T 2 A#) ~ 
1/(KA9) (where A9 is a typical angular scale for a change in a Lorentz factor. Thus, unless 
A# ~ 1/r 2 ~ 1/K, the term (delnR) 2 ~ 1/K 2 is of much higher order in 1/K and can be 
neglected. This express the fact the the lateral dynamics of strongly relativistic shock waves (in 
fact of any strongly relativistic motion) is "frozen out" by kinematic effects. 

If the shock is not strongly collimated, A# ^> 1/r 2 , we can neglect the factor dg lni?, different 
parts of the shock will propagate radially with a different Lorentz factor given by the driver or the 
external density inhomogeneity: 



This approximation is similar to Laumbach & Probstein (1969) approximation, which assumes 



radial motion. Thus, in the strongly relativistic case both Kompaneets and Laumbach-Probstein 



become equivalent, consistent with the conclusion of Shapiro (1979). 

In the strongly relativistic regime, K = T 2 1, the lateral dynamics of the flow is frozen out, 
unless the flow is extremely strongly collimated with A8 < 1/r 2 . Keeping (<%lni?) arbitrary and 
expanding in 1/K eq. Q takes the form 

dR\ 2 1 ( 1 



dt j 1 k yi + (d e inR) 2 ' (l2i 
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3. Discussion 

In this paper we considered the lateral evolution of non-spherical relativistic outflows. Con- 
trary to the commonly assumed fast lateral expansion, we find that unless the shape of the shock 
is extremely narrow, with the opening angle of the order of 1/T 2 , the lateral evolution is effectively 



frozen out by the highly relativistic motion of the shock. Thus, we confirm the conclusion of Shapiro 



( 1979 ) that highly relativistic shock propagate nearly ballistically. Our conclusion is broadly con- 
sistent with the results of numerical simulations showing very slow (logarithmic, e.g., ) lateral 
evolution. In contrast, the calculations of the afterglow emitted spectra have to be reconsidered 
accordingly. 

I would like to thank Dale Frail, Hendrik van Eerten, Andrew MacFadyen and Eli Waxman 
for discussion. 
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Fig. 1. — Geometry of the flow. The central source located at the origin produces anisotropic wind 
with luminosity depending on the polar angle 9. At a polar angle 9 the shock is located at radius 
R(9), while the wind direction (radial direction) makes an angle a with the shock normal 



